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A natural atom placed into a cavity with time-dependent parameters can be parametrically ex¬ 
cited due to the interaction with the quantized photon mode. One of the channels of such a process 
is the dynamical Lamb effect, induced by a nonadiabatic modulation of atomic level Lamb shift. 
However, in experiments with natural atoms it is quite difficult to isolate this effect from other 
mechanisms of atom excitation. We point out that a transmission line cavity coupled with a super¬ 
conducting qubit (artificial macroscopic atom) provides a unique platform for the observation of the 
dynamical Lamb effect. A key idea is to exploit a dynamically tunable qubit-resonator coupling, 
which was implemented quite recently. By varying nonadiabatically the coupling, it is possible to 
parametrically excite a qubit through a nonadiabatic modulation of the Lamb shift, even if the 
cavity was initially empty. A dynamics of such a coupled system is studied within the Rabi model 
with time-dependent coupling constant and beyond the rotating wave approximation. An efficient 
method to increase the effect through the periodic and nonadiabatic switching of a qubit-resonator 
coupling energy is proposed. 

PACS numbers: 42.50.Ct, 42.50.Dv, 85.25.Am 


I. INTRODUCTION 

Superconducting circuits with Josephson junctions can be used for quantum computing as qubits . These 
systems are macroscopically large, but they demonstrate quantum behaviour, which allows them to be treated as 
artificial atoms. Superconducting qubits can be integrated with microwave waveguides, while the photon field in these 
waveguides is quantized. Such systems give rise to the effects of Rabi oscillations [MQ], quantum feedback m, GHz 
photons emission [12] , while many-qubit systems form subwavelength quantum metamaterials [T3UT6] . A well known 
Rabi model [IT] is applicable for the description of qubit-photon quantum system [T]. 

Moreover, superconducting circuits integrated with microwave resonators provide a unique platform for an obser¬ 
vation of cavity quantum electrodynamical (QED) nonstationary phenomena which can hardly be studied in more 
traditional experiments. This happens because (i) resonator frequencies in superconducting systems are nearly six 
orders of magnitude lower than frequencies in optical systems, (ii) superconducting qubit-cavity systems demonstrate 
high and fast tunability of their main parameters. 

One of the examples of such nonstationary QED phenomena is the dynamical Casimir effect, predicted long time 
ago [18]. According to the initial idea, an ’empty’ space between two mirrors can emit photons due to vacuum 
fluctuations, provided these mirrors are rapidly moving with respect to each other. In order to observe this effect, 
one has to move mirrors with the speed approaching the speed of light, which seems to be unrealistic for a direct 
experimental realization. Eor this reason, there were several other proposals how the dynamical Casimir effect can 
be observed, see, e.g.. Refs. [IMB- In all these schemes, it was suggested that an effective boundary condition at 
’mirrors’ could be varied nonadiabatically, instead of moving massive ’mirrors’ themselves. It is remarkable that it 
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was a superconducting system, which led to the first observation of the dynamical Casimir effect [22l [23] . This was 
achieved by using an additional SQUID at the end of the waveguide in order to tune an effective boundary condition 
by varying in time the magnetic flux through the SQUID. 

A presence of an additional atom in a cavity with time-dependent parameters leads to other nonstationary QED 
phenomena, since an atom and a photon field can interact with each other [24ll28] . For instance, an atom can be 
parametrically excited, even if the cavity was initially empty. It was shown in Ref. [27] that there are several channels 
for such a process. In the case of a nonadiabatical modulation of cavity parameters, there are two channels of such 
a process. The most obvious mechanism is due to the absorbtion of Casimir photons. Another channel is due to the 
nonadiabatic modulation of atomic Lamb shift: virtual states of atom-photon coupled system are transformed into 
real states. This phenomenon can be called a dynamical Lamb effect. Note that the existence of a static Lamb shift 
for superconducting qubit energy levels was verified experimentally in Ref. [29]. In contrast to natural atoms, the 
effect can be significantly enhanced, since a strong coupling regime between the artificial macroscopic atom and the 
resonator is achievable. 

In practice, it is quite difficult to isolate the mechanism of atom excitation due to the dynamical Lamb effect from 
another channel due to the the absorbtion of Casimir photons. In order to overcome this difficulty, it was suggested 
in Ref. [27] that a single natural atom can be passed through the resonator consisting of two cylindrical cameras 
of different diameters which are characterized by two unequal Lamb shifts. When leaving one camera and entering 
another one, the dynamical Lamb effect can occur, so that an atom can be parametrically excited. A specific feature 
of this scheme is that, in this case, no Casimir photon appears, so that the dynamical Lamb effect can be, in principle, 
isolated. Unfortunately, a direct experimental implementation of this idea is rather difficult. 

The main goal of this paper is to suggest a superconducting qubit-resonator system as a suitable platform for 
the observation of the dynamical Lamb effect instead of a natural atom-cavity system. An astonishing property of 
superconducting qubits is a high tunability of their parameters in situ during an experiment [301 EH- For instance, 
one can modulate qubit excitation energy to manipulate qubits by using this degree of freedom [32]. Moreover, due 
to a recent technological progress, it becomes realistic to vary dynamically not only the resonator frequency, but 
also a coupling between the artificial atom (transmon qubit) and photon field with the frequency of a resonator or 
even faster. This quantity already can be tuned in sufficiently wide ranges by an auxiliary SQUID, see, e.g.. Ref. 
[33] and references therein. Using this additional degree of freedom, one can mimic a single atom passing from one 
camera to another one by a single switching of a qubit between the two resonators. Namely, this scheme can be 
mapped on a superconducting qubit at rest integrated with the resonator A, while at certain moment it is switched 
nonadiabatically to the initially empty resonator B. It is expected that this switching can be accompanied by the 
parametric excitation of a qubit. In this scheme, a channel of a qubit excitation through the absorbtion of real Casimir 
photons is suppressed, because in the weak coupling regime the perturbation of the cavity eigenmodes is negligible. 
Thus, instead of relying on a physical motion of an atom, one can simply tune a magnetic flux through a special 
auxiliary SQUID. 

Actually, a high tunability of superconducting circuit systems allows even for a significant simplification of this 
mapping. Instead of utilizing two resonators, one can use a single one. Indeed, one can consider a qubit initially 
uncoupled from the resonator, while at certain moment this coupling is nonadiabatically switched on leading again 
to a possibility of a qubit excitation and photon generation. Furthermore, a superconducting qubit-resonator system 
can be made in a strong coupling regime which must increase dramatically the probability of qubit excitation due to 
the dynamical Lamb effect. Note that an important condition to observe the dynamical Lamb effect is also that all 
switchings have to be performed nonadiabatically. This implies that the typical switching time to be smaller than the 
resonator frequency. Such a regime becomes possible due to the very recent progress [33] . 

In order to describe theoretically the dynamical Lamb effect in the system with the tunable qubit-photon coupling, 
we use the Rabi model. To treat the dynamical Lamb effect, we have to go beyond the rotating wave approximation [34] 
(RWA) in order to take into account counter-rotating processes. These are the photon creation with the simultaneous 
qubit excitation and the opposite process. Such processes in superconducting qubit-resonator systems are not at all 
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illusive, since they show up in experiments with strong-coupled systems, see, e.g.. Ref. [35]. We take into account 
counter-rotating wave terms by constructing a perturbation theory around the stationary RWA solution. In addition, 
we solve a time-dependent Schrodinger equation numerically. We find that our simple analytical formula for the qubit 
excitation probability gives excellent results in the case, when the interaction constant is switched nonadiabatically 
only once. 

Our treatment does not only yield the description of the system’s dynamics but also allows us to suggest a special 
trick which can be used to enhance a dynamical Lamb effect. Namely, we show that the qubit excitation probability 
can be dramatically increased provided that the coupling constant is switched on and off periodically with the period 
twice larger than the resonator frequency (each switching is nonadiabatical). Previously, similar periodic drivings have 
been proposed to enhance other cavity QED effects in a context of superconducting circuits, such as the dynamical 
Casimir effect [22l [23] or Lamb shift [36]. For such a periodic driving, our analytical perturbative result yields 
qualitatively correct description of the system’s dynamics. We also examine the statistics of photon states, generated 
due to the coupling constant dynamics. We show that a significant squeezing of these states can be achieved. This 
result can be of practical importance. 

We wish to stress that the main focus of this paper is not at all to develop a new method for superconducting qubit 
manipulation, but to suggest a scheme for the experimental realization of the dynamical Lamb effect. Note also that 
from the viewpoint of a quantum computing, the dynamical Lamb effect can be either positive or parasitic depending 
on a particular situation. For instance, it can lead to the undesirable qubit excitation. Therefore, the understanding 
and control of the dynamical Lamb effect is of importance both from the viewpoint of a fundamental physics and 
possible applications. 

This paper is organized as follows. In Section II, we present our idea and outline a theoretical model used. In 
Section III, we find qubit excitation probability for different regimes of qubit-resonator coupling dynamics. In Section 
IV, we address correlation functions for both the photon field and qubit degrees of freedom upon coupling energy 
dynamics. We conclude in Section V. 


II. MODEL 

There are different types of superconducting qubits. All of them represent a superconducting ring containing several 
Josephson junctions, but they differ from each other by ratios of their charging energies and Josephson energies. In 
order to study experimentally the dynamical Lamb effect, we need a qubit which allows for the dynamical tunability 
between its own degrees of freedom and photon field with the frequency of a resonator or even faster. Actually, these 
are charge qubits (transmons) and flux qubits that seems to be more suitable for such a task [33] . 

The cavity is a coplanar waveguide where the photon mode frequency uj is large compared to the characteristic 
energy of interaction between the qubit and photon degrees of freedom (weak coupling). This limit supposes that 
the resonator is operated in the single mode regime. In practice, in order to enhance the dynamical Lamb effect, a 
strong coupling regime might be more appropriate. However, it is much more difficult to address this regime within 
the analytical treatment. We, therefore, focus on a weak coupling regime and show that even in this case a significant 
effect is possible. 

The inductive coupling between a qubit and resonator photons is achieved by an auxiliary SQUID operating at the 
same GHz frequencies [33]. This opportunity is of a crucial importance for an observation of the dynamical Lamb 
effect. The coupling energy can be changed either by a single switching or via a periodic modulation. The latter 
approach will be shown to be much more efficient. 

The qubit-resonator system can be described in terms of the Rabi model naiMi, which is widely used in quantum 
optics. The total Hamiltonian of this model takes into account photons at the single mode uj and the qubit with the 
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bare excitation frequency e with the coupling V between them 

H = cua^a + 2^^^ ^ ^ 

where and a are secondary quantized boson operators of the photon field, while Pauli matrices as = 2cr+cr_ — 1, 
cr+, cr_ act in the space of qubit states. The operator V reads as 

V = g(a + a^)(a- + cr+), (2) 


where (a-ha^) and cr_ + cr+, up to numerical factors, are nothing but the electric field and dipole moment operators, 
respectively, while ^ ^ cj is a coupling const ant. Under this assumption, a perturbation of the cavity eigenmodes is 
negligible, so that a channel of qubit excitation via the absorbtion of real Casimir photons is not relevant. 

A decoherence k and relaxation rate T are of importance for real qubits, available in modern experiments. Typically, 
these two quantities turn out to be larger or of the same order as g. However, nowadays there is a steady and rather 
fast progress in decreasing /^, T, so we will concentrate on the limit g ^ T. 

In the case of a stationary system, a perturbation theory in operator V can be developed rather easily. It was shown 
in Ref. [27] that V produces contributions to the eigenenergies, which, in leading order, are independent on a photon 
level number but depend only on atomic level number. These contributions are interpreted as a Lamb shift, since 
these are virtual photons which are responsible for the effect. The approach of Ref. m is justified in the off-resonant 
regime |A| ^ where A = e — cj is a detuning frequency. It predicts the following shift for the atomic ground state 
energy: 


= - 


uj e 


(3) 


Actually, there is another approach to the problem, which has some advantages. Namely, the operator V can be 
split into the sum of two terms. 


V = V + V 2 , (4) 

where 

Vi = g(aa++ a^a-) (5) 

is the rotating wave contribution, conserving the excitation number, and 

V2 = g(a'^a++ acr-) (6) 

is the counter-rotating term, which breaks this symmetry, but conserves a parity. 

We will focus on the situation, when the qubit-cavity system is near the resonance, e uj, since in this case the 
dynamical Lamb effect is strongest. This resonant regime can be achieved by tuning qubit frequency e in situ via an 
external magnetic flux, threading the qubit loop. Near the resonance and if g is also independent on time, the fast 
oscillating term V 2 is usually omitted, while the resulting Hamiltonian is known to be exactly integrable. This well 
known solution is outlined in Appendix A. We have to take V 2 into account since it is precisely this operator, which 
is responsible for the stationary Lamb shift and, hence, for the dynamical Lamb effect as well. In the stationary case, 
V 2 can be treated as a perturbation. This approach is also utilized in Appendix A for the evaluation of the stationary 
Lamb shift in a resonant regime. Moreover, it will be demonstrated that a similar splitting of the Hamiltonian gives 
excellent results in certain nonstationary situations as well. 
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III. DYNAMICAL LAMB EFFECT 
A. Qualitative picture 

This section deals with the main subject of our paper. Namely, we consider a coupled qubit-resonator system, when 
the coupling energy g is time dependent. We assume that in the initial moment, the qubit and the resonator were 
not coupled, so that the initial state of the whole system was |0 |). When the coupling is switched on, higher-order 
photon states start to play their role, as well as a qubit excited state. In the case of a full resonance, this process 
can be illustrated by a scheme, presented in Fig. 1 . Apart from the bare ground state |0 |), there are two sets 
of doubly-degenerate excited bare states. These states feel each other through the operator V^ which is a sum of 
two contributions. The first one, Vi, conserves the excitation number, i.e., a total number of photons plus a qubit 
excitation number (0 or 1). The second one, ¥ 2 ^ changes the excitation number by 2 by adding/removing one photon 
and simultaneously switching qubit to the excited/ground state. 

Thus, we see that the bare ground state can be perturbed by V 2 only, which leads to the occupation of the |1 t) 
state with nonzero probability. In its turn, this latter state is connected by V 2 only with |0 1) and by Vi only with 
|2 I). By continuing this process, we see that (i) only the states with even excitation number are populated, (ii) the 
states with the same excitation number are linked only via Vi, (iii) the states with the excitation numbers differing 
by 2 are linked by V 2 only. 

Let us assume that tunable coupling g is not oscillating too fast. In this case, V 2 must be responsible to fast 
oscillations, while Vi corresponds to slow oscillations with frequencies essentially controlled by the Rabi frequency. 
Therefore, we may expect that V 2 can be treated perturbatively, although both Vi and V 2 are formally proportional 
to the same parameter g (a small parameter responsible for this feature in a resonant case is g/uj). This implies 
that the qubit excitation probability is going to experience oscillations with a frequency of the order of the Rabi 
frequency. Such oscillations are induced by processes between the states of the same excitation number. By solving 
time-dependent Schrodinger equation numerically, we will also show that this qualitative picture is indeed correct. It 
is also of importance to note that the photon-qubit dynamical coupling results in the occupation of all photon states, 
i.e., with both odd and even photon numbers. This feature, which can be probed in experiments, is in contrast with 
the dynamical Casimir effect, which leads to the occupation of the states with only even number of photons EH- 

Notice that slow Rabi-like oscillations show up even in the first order of the perturbation theory around the RWA 
solution. In contrast, when applying a perturbation theory around bare states, all terms of an infinite expansion in 
V must be taken into account in order to recover such oscillations. This is one of the advantages of the theoretical 
approach utilized in the present paper. 


B. Time-dependent perturbation theory 

As a zero-order approximation, we take the stationary RWA Hamiltonian, in which the coupling constant is equal 
to the time average value g = {g{t))t (averaging is performed after the coupling is switched on) and remaining 
nonstationary terms are treated as a perturbation. This perturbation contains counter rotating wave operator V 2 
and a rotating wave contribution, which is due to deviations of g from its time averaged value. Namely, the full 
Hamiltonian is represented as 


H{t) = HiiWA+V{t), 


( 7 ) 


where 


Hrwa = + -e(l + 0 - 3 ) + g{aa+ + aV_), 

V{t) = g{t){a^a+ + acr_) + {g{t) - g){aa+ +aV_). 


( 8 ) 
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FIG. 1: A scheme illustrating an internal structure of the full wave function projected into states of the qubit-waveguide 
system without coupling (see in the text). 


Let us now consider the evolution of the wave function = 0) = |0 |) assuming that a coupling is switched on at 
t = 0. Then, at the moment t > 0, the wave function is given by the standard expression 

tp{t) = I -i j 

V 0 

where T stands for the time ordering operation. In our calculations we take into account only first order terms in 
T-exponent series expansion in Eq. This approximation is further justified by a direct comparison with the result 
of the solution of the time-dependent Schrodinger equation. 



C. Driving by a single switching of a coupling constant 

We now study the evolution of the wave function in a full resonance regime, e = cj, after the sudden switch of the 
coupling which then stays constant in time. Thus, we have g{t) = g0{t) and g = g, so that the second contribution 
in operator E(t), as given by Eq. (|^, vanishes. Using Eq. , we obtain 

ip{t) = lOJ,) - — (e“ 2 “*cos\/ 25 i- 1)|11) - smV2gt\2 i). (10) 

2 c(j 2 c(j 

In this approximation, the qubit excitation probability as a result of the instantaneous switching, is given by 

2 

w^{t) = ^3 + cos 2V2gt — 4 cos 2 co’t cos \/ 2 ^t^ . (11) 

This probability corresponds to the dynamical Lamb effect and includes both slow and fast oscillations. In the 
superconducting systems with weak coupling between photon field and qubit, the probability is as small as g‘^luj‘^. 
Note that the amplitude of this quantity is proportional to the square of the change of the static Lamb shift, in 
accordance with Ref. m- 

The number of generated photons as a function of time nph{t) = a\ 2 p{t)) can be also readily extracted from 

Eq. ( [Tq| ). We would like to stress that these photons appear as a result of the dynamical Lamb effect. Our results for 
w^{t) and nph{t) ai uu/g = 20 are plotted in Eig. 2 for 0 < t < 1.2 Tr, where the Rabi period is defined as Tr = ir/g. 
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In order to justify these results, we also solved Schrodinger equation numerically in the weak coupling regime. The 
basis was truncated by typically 10 photon states. An accuracy was verified by increasing the number of such states 
and comparison with results for a smaller basis. The results of such numerical computations show excellent agreement 
with our perturbative results for both w^{t) and nph{t). The relative discrepancy is essentially within few percent. 
Such an excellent agreement justifies our qualitative understanding of the system dynamics, discussed above. 


w-rit) 




FIG. 2: The qubit excitation probability (o-) 9^nd the photon number nph{t) (b) as functions of time after a single 

switching of the qubit-waveguide coupling constant at cu/g = 20. 


D. Parametric driving 

Single switching of the qubit-resonator coupling results in the low probability for the qubit excitation. In order to 
enhance it, one can use more sophisticated modulations of ^ as a function of time. We, however, should keep in mind 
that the dynamical Lamb effect can be isolated only in the case, when this coupling changes nonadiabatically, so that 
the only one way is to switch it off and on instantaneously. 

The excitation probability can be dramatically increased at the regime of a parametric resonance, when the driving 
is induced by high frequency periodic switching of the qubit-resonator coupling. In particular, we suggest applying 
2co’-periodical switching on and off, namely, g{t) = g0{cos2ujt). Hence, we take g = g/2 in the Hamiltonian splitting, 
as given by Eq. (U). In order to evaluate qubit excitation probability, we again limit ourselves to the first order 
approximation in Eq. 












































































We again assume that uj ^ g. Therefore, fast oscillating terms of the form can be replaced by their average 

values when performing integration in Eq. ([^, while corrections to this approximation are negligible as 

g/cj. Under this replacement, we find 


'0(t) |0 I) H - e 

TT 


The corresponding excitation probability is 


— 2iujt 


1 - cos |2 4 ,) - zsin ^ | 11 ) 


w^{t) ^{1 — cos V2gt). 


( 12 ) 


(13) 


This result describes an approximate Rabi-like oscillatory behavior of the qubit excitation probability via the dy¬ 
namical Lamb effect in the regime of a parametric resonance. Note that frequences of oscillations of qubit excitation 


probability and mean photon number differ by a factor of 2 , according to Eq. ( 12 ). Regardless of the relation between 
amplitude of the switching g and cj, the maximum probability reaches the finite value 4 / 7 r^ on a characteristic time 
scale ^ Xjg^ which depends on the amplitude of the inverse coupling strength only. In contrast to the probability 


in the regime of a single switching the time dependence of w^{t) in parametric resonance regime does not 
contain high-frequency oscillations at uj. Moreover, the maximum probability is not small anymore. This last feature 
is of great importance for the possibility of an experimental investigation of the dynamical Lamb effect. 

We also wish to stress that the precise value of the maximum excitation probability, as given by Eq. ( pT] ), is rather 
relative. Namely, if we take into account the next term of the perturbative expansion in Eq. ([^, the prefactor 2 / 7 r^ is 
going to change, staying nevertheless independent on In addition, an admixture of Rabi-like harmonic with higher 
frequency will appear, which corresponds to oscillations between the two states with the number of excitations equal 


to 4. Nevetheless, Eq. (13) yields correct qualitative description for as further evidenced by a comparison with 

fully numerical results. A lack of a quantitative agreement is due to the fact that our picture based on the separation 
of fast and slow degrees of freedom is not so accurate anymore because of a high-frequency modulation of g. 

Next, we solve the time dependent Schrodinger equation numerically. We indeed see Rabi-like oscillations for both 
the qubit excitation probability and the mean photon number, as shown by Eig. 3. Erequencies of these oscillations 


are correctly described by Eq. (13). The estimated period of w^{t) oscillations corresponds to the time scale V^jg 
appearing in the RWA Hamiltonian in the two-excitation solution. We also found that maximum w^{t) and photon 
number nph{t) are only slightly sensitive to the amplitude of the coupling strength g in the weak interaction regime, 
which is again in a qualitative agreement with the simple analytical treatment. At the same time, the value of 


maximum as predicted by Eq. (13), is nearly twice underestimated compared to the exact (numerical) solution. 

Thus, both w^{t) and nph{t) demonstrate a universal behavior. Namely, by rescaling the time, they can be cast in a 


universal form within a few percent accuracy. This feature, which is predicted qualitatively by Eq. (13), is in a full 
agreement with the results of our numerical computations. 

Recently, it was suggested theoretically that the Lamb shift of a qubit can be greatly enhanced via periodical 
driving of a qubit by the external classical electromagnetic field [36]. This idea has some similarities with our idea 
of parametric pumping. However, in contrast to Ref. [36], we do not need an external classical field which directly 
affects a qubit and modifies its Lamb shift. Instead we modulate only the qubit-resonator coupling, while qubit is 
excited and photons appear as a result of the interaction between a qubit and resonator. 

Nevertheless, both ideas indicate that counter-rotating processes in a weak coupling regime can be significantly 
enhanced by using various types of a periodic driving. The same strategy was used to increase a dynamical Casimir 
effect which finally has led to its experimental observation [22j|23]. Thus, periodic driving, in general, can be exploited 
to enhance various nonstationary cavity QED effects. 
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FIG. 3: The qubit excitation probability w^{t) (a) and the photon number nph{t) (b) as functions of time after the beginning 
of parametric driving at = 20. 


IV. CORRELATION FUNCTIONS 


In addition to the qubit degrees of freedom, our system has also photon ones. They can be characterized by 
various correlation functions. For instance, one can use squeezing parameters defined in a standard way as Ap = 
I a/(( a — at)^) — {a — at)^ and Ax = ^ {{a + ai)^) — (a + ai)^. For the stationary system, these two parameters can 

be evaluated approximately using a perturbation theory around RWA solution (see Appendix A). The wave function 
yields the following squeezing parameters (in the first order by g) 


Ap = Ax = 


1 

2 



294-^ 


(14) 


It is easy to see that Ap, Ax > 1/2, since < 0, which implies that this state is not squeezed. 

Let us now consider squeezing in the case of a parametric modulation of the qubit-resonator coupling constant. 
An accuracy of the analytical solution in the first order of perturbation theory is not sufficient for the quantitative 
analysis. We therefore restrict ourselves to the numerical solution of the time dependent Schrodinger equation taking 
into account nearly 30 photon states. Numerical solution shows the realization of the squeezed state under parametric 
driving of a particular form g(t) = g0{cos{2ujt)). We observe non-periodic oscillations of the squeezing parameter Ax. 
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Ax{t) 



FIG. 4: Numerical result for the squeezing parameter Ax(t) at oj/g = 20 within 20 Rabi periods after the beginning of the 
parametric driving. Only lower and upper envelope curves are shown, while Ax{t) experiences fast oscillations between them. 


Figure 4 shows the lower and upper envelope curves for the evolution of Ax within first 20 Rabi periods after the 
switching on the driving g{t) for uj/g = 20. These envelope curves experience slow oscillations, while fast oscillations 
of Ax occur between them. The fast oscillations are not shown in Fig. 4. 

The results again reveal the universality of the envelope curves oscillations, regardless of the ratio uj/g. The 
minimum value of Ax at particular moments of time, according to Fig. 4, is ~ 1/4, which indicates a sufficiently 
strong squeezing of a generated photon state. By comparison Fig. 4 and Fig. 3, we see that oscillations of Ax are 
correlated with oscillations of nph{t). 
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FIG. 5: Numerical result for photon correlation functions at to = 0.5 Tr (a), 4Tr (b), STr (c) at uj/g = 20 within 

20 Rabi periods after the starting of the parametric driving. 


The first-order and second-order correlation functions can be expressed through the field operators in a Heisenberg 
picture as 

G (toTj - (at(to)a(to)>’ 

_ {aHto)aito)aHt)ait)) r\ 

V / {a^to)a{to)a^to)alto)) ’ v / 

These two quantities are calculated using a numerical solution of the time-dependent Schrodinger equation. Typical 
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FIG. 6: Numerical result for photon correlation functions at to = 0.5 Tr (a), 4Tr (b), and STr (c) ai u/g = 20 

within 20 Rabi periods after the starting of the parametric driving. 


results are presented in Figs. 5 and 6. Both correlation functions demonstrate non-periodic oscillations with slow 
decay as functions of t — to, where to corresponds to the starting point and the correlation is probed at t > to- By 
comparing these two figures with Fig 3, it is easy to see that the correlation functions at small to generally follow the 
field revivals as a function of t, while oscillations are not so pronounced at larger to- The existence of slow oscillations 
is clearly linked to the interplay between the qubit and photon field, which is accompanied by the qubit excitation 
due to the dynamical Lamb effect. 

We also introduce first-order and second-order correlation functions for qubit degrees of freedom (qubit excitation 
number) as 


Gq\to,t) 


(cr+(to)o-_(t)) 

(cr+(to)cr-(to)) ’ 


^(2)/. _ (cr+(to)cr_(to)o-+(t)cr_(t)) 

<7 V O 5 ; (cr+(to)cr_(to)cr+(to)cr_(to)) ’ 


(16) 


which are again calculated numerically. The results are shown in Figs. 7 and 8 (for the same set of parameters as for 
Figs. 5 and 6). The behavior of these correlation functions generally resembles the behavior of similar quantities for 
the photon field. The dynamics of and follows the dynamics of 

However, the amplitude of oscillations of G^q^ is rather small, while this quantity does not experience a decay as a 
function of t — to. This is related to the two-level character of qubit energy spectrum which has to be contrasted with 
the unbounded photon spectrum. Moreover, for the considered range of parameters, photons occupy several levels, 
while qubit mostly stays in its ground state. This latter feature leads to only small oscillations of G^q\ Nevertheless, 
slow oscillations of correlation functions for the photon field, which can be probed in experiments, are definitely related 
to slow oscillations of correlation functions for qubit degrees of freedom, as well as to the qubit excitation probability. 

Thus, slow Rabi-like oscillations of various characteristics of a photon field can be used in experiments as a signature 
of the dynamical Lamb effect in the systems with tunable photon-qubit coupling. Indeed, if the resonator was initially 
empty and qubit was uncoupled from it, no such oscillations can exist, of course. If a coupling energy is modulated, 
qubit can be excited leading to photon creation via counter-rotating processes, which can be seen in experiments. In 
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FIG. 7: Numerical result for qubit correlation functions Gq^\to,t) at to = 0.5 Tr (a), 4 Tr (b), STr (c) ai oj/g = 20 within 20 
Rabi periods after the starting of the parametric driving. 


contrast to the dynamical Casimir effect, photons occupy not only even states. 


V. CONCLUSIONS 

We suggested that a coupled superconducting qubit-resonator system can be used for the observation of the dy¬ 
namical Lamb effect. This effect was initially introduced for the natural atom in a varying cavity m- It leads to 
the probability for the atom to be excited solely due to the nonadiabatic modulation of atomic levels Lamb shift. 
Although this effect was initially predicted for natural atoms, it is quite problematic to observe it in experiments 
with such atoms, since this effect must be carefully separated from other mechanisms of atom excitation. To this 
end, it was proposed in Ref. [27] to pass a single atom through the resonator consisting of two cameras with different 
diameters. When passing the cameras, an atom can be excited and a photon can be generated. We suggest a non¬ 
trivial link between this scheme and a system with a single superconducting qubit (macroscopic artificial atom) at 
rest integrated with the resonator in such a way that the coupling can be tuned dynamically. By switching on this 
coupling nonadiabatically, one can excite a qubit and generate a photon, even if the resonator was initially unoccupied 
by photons. Such a dynamical coupling between the qubit (transmon) and the resonator is already achievable by 
utilizing an auxiliary SQUID [33]. Thus, we suggest that this scheme can be used for the experimental realization of 
the dynamical Lamb effect. 

We presented a theoretical model for the description of the system dynamics upon a modulation of the coupling 
energy and for the predictions of a method to increase the effect. Our treatment is based on the dynamical Rabi 
model in which we take into account counter-rotating terms, since they are responsible for the dynamical Lamb effect. 
The dynamical problem was solved both analytically using a perturbation theory around the RWA stationary solution 
and numerically by integrating time dependent Schrodinger equation. 

We found that a single switching of the coupling constant results only in a low probability for the qubit to be 
excited when considering weak coupling regime, as well as a full resonance between the qubit and the cavity single 
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FIG. 8: Numerical result for qubit correlation functions at to = 0.5 Tr (a), 4 Tr (b), and 8 Tr (c) aX uj/g = 20 within 

20 Rabi periods after the starting of the parametric driving. 


mode. After the switching on, this probability experiences slow Rabi-like oscillations. We then suggested a much 
more efficient method to enhance the effect which is a periodic nonadiabatical switching of the coupling energy. We 
also studied a statistics of photon states generated upon the coupling constant modulation. We found that these 
states can experience a strong squeezing. 


Acknowledgments 


The authors acknowledge useful comments by A. V. Ustinov, A. L. Rakhmanov, V. P. Yakovlev, and L. V. Bork. 
This work was supported by RFBR (project no. 15-02-02128). W. V. P. acknowledges the support from the Russian 
Science Support Foundation. Yu. E. L. is supported by Program of Basic Research of HSE. 


Appendix A: Perturbation theory around RWA solution 


The stationary Rabi Hamiltonian can be split as 

H = Hrwa + V2. 

The first part Hrwa is the Hamiltonian in the rotating wave approximation. 


Hrwa — ida^a + x^(l + crs) + Vi 


1 


(Al) 


(A2) 


which is known to be exactly integrable. 

The set of the ” dressed” eigenstates and corresponding eigenenergies of Hrwa is given by the following expressions 


i’n'’ = ai '>\ni)+l3^ V - 11), 

Xn ^ k i) + k -11) , e; 


XO) 

n,'ip 

( 0 ) 


(A3) 








14 


where \n \) and |n — 1 t) form a set of ’’bare” states of the decoupled qubit and photons. Coefficients in (A3) read 
as (n > 1 ): 


a 


(±) 


gVn 


+ 2g^n ± AJ ^ + g^n 


(A4) 




f±J^+g^n 


^+2g^n-AJ^+g^n 


A bare ground state |0 |) is not affected by Vi and remains the same as at ^ = 0. Thus, there is only one state 
corresponding to n = 0 , which is 


4°^ = 104 ,). 


The eigenenergies of Hrwa are 


A /A2 

=ujn + — - \ + g^n, n > 0 , 


(AS) 


^ri,x 


A /A2 

+ — + Y — + gM, n>l. 


Using standard perturbation theory, we find that the wave functions within the first order contribution in V 2 read 
(n > 1 ) 


^„/(0) , /^n+2^^ ^ ^/n 1 fQ) ^ a/tI 1 ('q) ^ VTl + 1 ('q) 


4^^ = V’r +51 - C-(2 + 'r'_ • -Xn-2 + ■ _ , 

\ ^n—2,'ip ^n,'ip ^n+2,4 ^n—2,x ^n+2,x 


Xn +2 ^ • 

(A 6 ) 


( 1 ) ( 0 ) I f *^n-2A \fn T ( 0 ) 4+2®” 4^ + 1 ( 0 ) 1 ,( 0 ) 4+4” 4^ + 1 , ( 0 ) 

XV = XV + 51 ^—x „-2 + —x „+2 + ———- 4-2 + , - 4+2 


\ ^n,x 2 ,x ^'^+ 2 ,x 

For the modified ground state (n = 0), we obtain 


^n,x ^n—2,4 


^n,x ^n+2,4 


41 ) = |0 ;) - ^12 I) + 4-4 t))■ 

^2,4 


(A7) 


(AS) 


Using the above strategy, one can also find corrections in V 2 to the eigenstates of Hrwa- For instance, a first 
non-zero contribution to the ground state energy is determined by the second order correction and is given by the 
negative value 


(2) _ g^l4~4 ^ -V(l-A/4A2 + 8g2) ^ 

^2,x 3cj + e — y^A^ + 8g‘^ 


^ 0,4 ~ 


(A9) 


This is the Lamb shift value, where all the processes, generated by the rotating wave term Vi in the Hamiltonian, are 
included in zero-order approximation. This scheme of the expansion with respect to the operator V 2 around Hrwa 
regularizes the perturbation theory against 1 /A divergencies. It is applicable in the opposite resonant case |A| ^ g^ 
where the perturbative result of Ref. [27] is not adequate due to divergent corrections ^ g/A to the wave functions 
and eigen energies. The divergencies originate from the degeneracy of the self-mode and qubit energy. 
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It is of importance to note that Eq. (A9) does not match the off-resonant result (§ that means that this V 2 - 
expansion is most reliable in the resonant limit A <C 5 ^, despite that the detuning A is formally not restricted in this 
approach. In order to recover the off-resonant result we have to consider contributions from higher orders of V 2 . 
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